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Effects of sliort-range correlations on the Coulomb screening, the phonons, and the pairing 
interactions are examined in electron-phonon systems. First, we derive a model Hamiltonian 
of Coulomb interactions which includes both the long-range part Vq and the short-range part 
U. It is found from the expression of the dielectric function that the strong on-site correlations 
weaken the Coulomb screening. Secondly, we examine the screened phonons and the interaction 
mediated by phonons. In a consistent picture, we derive an expression of the effective interac- 
tion which includes (1) the screened Coulomb interactions, (2) the pairing interactions mediated 
by phonons, and (3) the effective interactions mediated by spin and charge fluctuations. It is 
rewritten in a form of a summation of (a) the effective interactions of the pure Hubbard model 
without the long-range Coulomb interactions, and (b) the phonon-mediated interactions plus 
screened Coulomb interactions with corrections due to both U and Vq. Thirdly, we derive an 
effective Hamiltonian analogous to the BCS Hamiltonian. Fourthly, for some typical values of 
parameters, we obtain the ground state phase diagrams. It is found that spin-triplet supercon- 
ductivity mediated by phonons occurs when the short-range electron correlations are sufficiently 
strong, and the Coulomb screening is sufficiently weak. We estimate the orders of the transition 
temperatures when the triplet superconductivity occurs. The obtained values are realistic for 
existing candidates of the triplet superconductors as the order of the magnitudes. The possible 
relevance of the phonon-mediated interactions to the heavy fermion superconductor UPts and 
the layered superconductors such (TMTSF)2X and Sr2Ru04 are briefly discussed. 

KEYWORDS: phonon-mediated pairing interactions, Coulomb screening, electron-phonon system, over- 
screening, spin and charge fluctuations, exotic superconductors, spin-triplet pairing, heavy 
fermion superconductors, organic superconductors, ruthenate superconductors. 



1. Introduction 

It is widely known that overscreening of the repulsive 
Coulomb interactions between electrons by ions gives rise 
to attractive interactions for low frequencies |a;| < uju-^ 
This dynamical effect is considered to be the origin of the 
pairing interactions mediated by phonons, which induce 
the superconductivity in metals. At the same time, the 
long-range Coulomb interactions are screened into short- 
range interactions, and consequently the energies of the 
optical phonons of short wavelengths are reduced into 
the accoustic type ujq ~ |qf|. 

On the other hand, the screening of electronic and ionic 
charges influences the momentum dependence of the at- 
tractive interactions through the dielectric function. If 
there were not the screening, the coupling constant be- 
tween electrons and optical phonons diverges as |q| ^ 
due to the long-range nature of the Coulomb interactions. 
However, in actuality, such divergence is suppressed by 
the screening, but the coupling constant remains momen- 
tum dependent. 

Owing to such momentum dependences of the 
electron-phonon coupling constant and the phonon en- 
ergy, the resultant pairing interaction mediated by 



phonons depends on the momentum transfer. Hence, it 
includes the anisotropic components, when it is expanded 
in terms of appropriate basis functions. Usually, such 
anisotropic components are ignored, because they are 
considered to be much smaller than the isotropic com- 
ponent in most cases. In fact, the superconductivity in 
the ordinary metals, which are probably mediated by 
phonons, are considered to be of s-wave pairing. 

However, as Foulkes and Gyorffy pointed out,^ when 
the short-range Coulomb interactions are strong enough 
to suppress isotropic pairing, anisotropic pairing can 
be induced by the subdominant anisotropic interaction. 
They have examined possibility of spin-triplet supercon- 
ductivity by this mechanism in metals such as Rh, W, 
and Pd. 

Abrikosov discussed high-Tc cuprates, and proposed 
that such momentum dependence of the phonon- 
mediated pairing interaction could give rise to ex- 
tended s-wave superconductivity,'^ although it seems 
from some experiments that the cuprates are d-wave su- 
perconductors. Following his theory, Bouvier and Bok 
calculated the superconducting gap function and ob- 
tained anisotropic momentum dependence.'* Friedel and 
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Kohmoto, and Chang, Friedel and Kohmoto examined 
the d-wave pairing interactions mediated by phonons in 
the cuprates.^'^ The present author and Kohmoto have 
proposed^ that in the presence of coexisting ferromag- 
netic long-range order, a subdominant triplet state oc- 
curs after the singlet states are suppressed by the strong 
exchange field. We have discussed a possibility of this 
mechanism in UGe2.*'^ 

It is easily verified that the momentum dependence 
of the interaction is stronger for weaker screening. For 
example, in the Thomas-Fermi approximation, the cou- 
pling constant has a peak with the width of the inverse 
screening length in the momentum space. Therefore, in 
the interactions, the ratios of the anisotropic compo- 
nents to the isotropic component increase^" when the 
screening becomes weaker. There are some possible fac- 
tors which weaken the screening effect, such as crystal 
structures, low density of charge carriers, and low den- 
sity of states. The present author and Kohmoto have ex- 
amined layer structures as the weakening mechanism, 
because some candidates of the triplet superconductors 
have layer structures as wo shall discuss below. It was 
found that when the layer interval increases, the screen- 
ing becomes weaker, and subdominant triplet pairing 
interaction could more easily overcome the dominant 
s-wave pairing interaction with an asist of the on-site 
Coulomb repulsion. 

It should be noted^'^^ that the on-site Coulomb in- 
teraction is specific to the electron systems in solid, and 
does not appear in the electron gas model in the contin- 
uum space. The magnitude of the on-site Coulomb en- 
ergy U depends on the profile of the Wannier functions. 
Since such an on-site interaction is constant in the crys- 
tal momentum space, it only suppresses isotropic pair- 
ing, and consequently it favors subdominant anisotropic 
pairing. 

With respect to application, we are motivated by 
recently discovered exotic superconductors, such as 
Sr2Ru04, organic superconductors, and heavy fermion 
superconductors. In the Sr2Ru04 superconductor, spin- 
triplet pairing has been suggested by a Knight shift mea- 
surement^^ and a iiSR experiment. Many theories have 
been proposed on this material, and there arc some con- 
troversies especially on the momentum dependence of the 
superconducting gap function. Particularly, some non- 
phonon mechanisms have been considered.^'* but there 
has not been any evidence at present. On the other hand, 
the Sr2Ru04 superconductor exhibited the isotope effect 
coefficients with unusual dependence on the impurity (or 
oxygen defect) concentration and a reverse isotope effect 
in clean samples. The formar can be explained by as- 
suming internal transition of superconductivity,*® if the 
pho no n- mediated pairing interaction exists. For the lat- 
ter result, it has been shown^'' that within the weak cou- 
pling theory, the reverse isotope effect indicates that the 
Coulomb interaction between electrons, in total effect, 
works as repulsive interactions against the superconduc- 



tivity, unless the system is in the very vicinity of any 
magnetic instability. At present, any pure nonphonon 
theory could not explain these experimental facts of the 
isotope effect. 

Another candidate of the present mechanism of 
triplet superconductivity is the family of (TMTSF)2X 
compounds. From the experimental results at the 
early stage after the discovery of superconductivity in 
these compounds, including nuclear magnetic resonance 
(NMR)**' and the phase diagrams,^" c?-wave supercon- 
ductivity mediated by antiferromagnetic spin fluctuation 
has been studied in these compounds.^*' However, re- 
cent experimental studies suggest spin-triplet supercon- 
ductivity in (TMTSF)2PF623 and (TMTSF)2C104,2*'2^ 
although some of the experimental results have not been 
explained by triplet pairing.^® We will examine this prob- 
lem in a separate paper in details, applying the present 
mechanism of triplet pairing to these compounds. Only 
by the phonon mechanism, it might be difficult to re- 
produce the pressure dependence of the superconduct- 
ing transition temperature Tc observed in these com- 
pounds. The pressure dependence may be explained 
by taking into account the contribution from the spin 
fluctuations to the pairing interactions. The pairing in- 
teractions mediated by the spin fluctuations include the 
attractive triplet components. ^^'^^ 

We also consider the heavy fermion superconductors, 
such as UPta,^^'"^" as candidates in which the pairing in- 
teractions mediated by phonons may contribute to spin- 
triplet pairing. Since the width of the effective band is 
extremely narrow, s-wave state could not avoid the on- 
site Coulomb repulsion by the retardation effect. If one 
use the Thomas-Fermi approximation in the estimation 
of the screening length, it becomes much smaller than 
the lattice constant. However, it is not justified for a 
phenomena of such a small length scale. In actuality, the 
dipole field due to the charge distributions in the unit 
cell could not decay so rapidly within the scale of the 
lattice constant. Hence, the electron-ion interaction is 
not necessarily limited on the each lattice site. There- 
fore, it may contribute to the anisotropic pairing inter- 
actions. Whatever the pairing mechanism is, since the 
heavy mass renormalization significantly reduces the ef- 
fective coupling constant of the pairing interaction, the 
bare pairing interaction needs to be very large for the 
observed Tc to be reproduced. 

In general, it is reasonable that the phonon-mediated 
pairing interaction could contribute to superconductivity 
to some extent in most cases, even when the resultant gap 
function is anisotropic. In fact, in almost all superconduc- 
tors including exotic ones, non-vanishing isotope effects 
have been observed. It is reasonable to take them due 
to attractive contributions from phonons to the pairing 
interactions. It is likely that the momentum average like 
{j* (k)Vph{k — k')'-f{k')) is negative, since the phonon- 
mediated pairing interaction Vph{q) have a negative peak 



J. Phys. Soc. Jpn. 



Full Paper 



Hiroshi Shimahara 3 



around q = and |7(fe)| > 0, where 7(fc) denotes the 
function that expresses the momentum dependence of 
the gap function A(fc) oc j{k). Further, it is Ukely that 
the isotropic component is dominant, and the second and 
third-dominant ones are anisotropic and of odd and even 
parities in the momentum space, respectively. 

In particular, in the spin-triplet superconductors, the 
phonon-mediated pairing interaction can be the main 
mechanism of the anisotropic superconductivity, al- 
though another mechanism might assist it at the same 
time.^^ By an analogy to the superfluidity in liquid 
^He,^^ the paramagnon mechanism is sometimes consid- 
ered to be responsible for the spin-triplet superconduc- 
tivity. However, the situation is very different from the 
electron systems in the solids. The liquid ^He does not 
have lattice vibrations, and its transition temperature is 
of the order of ImK.^^ Therefore, we should be careful 
when we use the analogy. 

In the anisotropic singlet superconductors, the third- 
dominant component of the phonon-mediated interaction 
may contribute to the superconductivity to some extent. 
However, we presume that it could not be the only ori- 
gin of the pairing interaction. If the pairing interaction 
of nonphonon origin is negligible, the system should un- 
dergo a transition to triplet superconductivity at a higher 
temperature, since the third-dominant even-parity (spin- 
singlet) component is smaller than the second-dominant 
odd-parity (spin- triplet) component. 

In this paper, we examine the effects of the short- 
range (on-site) correlations on the Coulomb screening, 
phonons, and the pairing interactions, taking into ac- 
count the long-range Coulomb interactions at the same 
time. The main purpose of this paper is to devolop a 
general theory, although we are motivated by the su- 
perconductors discussed above. In particular, we derive 
a general effective Hamiltonian in a unified framework 
in which the short-range and long-range parts of the 
Coulomb interactions are treated consistently. The effec- 
tive Hamiltonian includes (1) the screened Coulomb in- 
teractions, (2) the phonon-mediated interactions by the 
screened electron-phonon interactions, and (3) the effec- 
tive interactions mediated by spin and charge fluctua- 
tions. In appropriate limits, it is reduced to the effective 
Hamiltonians examined by many authors so far.^"^'^^ In 
particular, in our previous papers, the on-site Coulomb 
energy U has been treated as a parameter independent 
of the phonon-mediated pairing interaction, but in prac- 
tice it also modifies the pairing interactions. We will also 
clarify it in this paper. 

In §2, we examine the model of the electron-phonon 
system. Starting from a basic model of the electron- 
phonon system, we derive a Hamiltonian which includes 
both the short-range part U and the long-range part Vq 
of the Coulomb interactions. In §3, we apply the random 
phase approximation (RPA) to the Hamiltonian obtained 
in §2, and examine the screening effects on the Coulomb 
interactions and phonons. The approximation is reduced 



to that in the electron gas model if we put U = and 
Vq ^ 0, while that in the pure Hubbard model if we put 
U ^ Q and Vq = 0. In §4, we examine the two-particle 
vertex part which contributes to superconductivity. We 
derive the general form which includes corrections due 
to U as mentioned above. In §5, we derive an eff'ective 
interaction within the weak coupling theory from the re- 
sults of §4. We examine anisotropic superconductivity on 
the basis of the effective interaction. Expressions of the 
transition temperature and the isotope effect coefficient 
are derived. In §6, we apply the effective model to some 
typical cases to clarify essential aspects of the present 
mechanism. The phase diagrams and the transition tem- 
peratures are obtained. The last section §7 is devoted to 
summary and discussion. 

2. A Model of the Electron-Phonon System 

In this section, we derive a model that we will exam- 
ine in this paper. We start with the general form of the 
coupled electron-phonon Hamiltonian 

H = Ho + ife-ph + He (2.1) 

with 

^^0 = ^ ekcl^Cka + ^ nqxbl^bqx (2.2) 

q\ 



ka 



kqaX 



(2.3) 



^^ = 1 Yl ^kX,k,kA^a<^k,aci^,,Ck,a'. (2.4) 



Here, Cfco- denotes the electron operator of the Bloch state 
with the crystal momentum k and the spin a in the rel- 
evant electron band, while bqx donotes the bare phonon 
operator of the phonon with the momentum q and the 
phonon-mode A. Since we consider the longitudinal op- 
tical phonons in this paper, we omit the phonon-mode 
suffix A and the summation over it for simplicity. Thus, 
we put flqx = Op, where ftp denotes the ionic plasma 
frequency. We consider the positive background charge 
of the atoms within the jellium model. ^ Therefore, we 
put the matrix element of the electron-phonon interac- 
tion Mkqx = Mq which satisfies 



2Ml _ 



(2.5) 



where VccW denotes the unit cell volume. The present 
model might not be very accurate, but since we need an 
order estimation at most for our purpose, we use this 
model in this paper for simplicity. 
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The coupling constant V'^jfe^fegfe^ is expressed as 



and 



in terms of the Bloch wave functions 



(2.6) 



(2.7) 



where Uk{r) satisfies the periodicity condition Uk{r + 
R) — Uk{r) for any lattice vector R, according to the 
Bloch's theorem. 

2.1 An expression in the tight-binding model 

Now, we rewrite the Coulomb interaction Hq in the 
tight-binding model. When the tight-binding approxima- 
tion is valid and the electron wave function has a large 
amplitude only near the atomic site, we may assume that 



R 



Hr-R), 



(2.8) 



where (/)(r — R) denotes a localized orbital wave func- 
tion near the atomic site at R. This approximation is 
equivalent to putting 



Uk{r) 



1 



(2.9) 



Since we may put r « i2 in eq. (2.9), it is equivalent also 
to Ufe(r) « Af-V2^^^(r-i2).If we define the Wannier 
function 

wir -R) = J^Y^e-'^-^Mr), (2.10) 

we obtain w{r) = (f){r) from eq. (2.8). 
Within this approximation, we obtain 



TACTCT 



/dVdV \^{r-R) 

RR' ^ 



X U(r'-i2')Pe"'^''^"'"=^-^T 



-i{k^—k^)-R' 



(2.11) 

where we have omitted terms including the factor of the 
form — Ri)(j)(r — R2) with Ri 7^ R2, since it is small 
from the assumption mentioned above. In the integrand 
of eq. (2.11), wo could replace r and r' in e^/jr — r'| 
with R and R' , respectively, due to the factors of \(p{i — 
i?)p|0(r' - H')P, except when R ^ R' . We could not 
make this replacement when R= R', because e^/jr — r'| 
diverges when 1 — > r'. Therefore, we obtain 



^fe'lfe2fe3fe4 = {U^aa' + Vki-k2)^ki-k2,k3-k4, 

with Sera' = 1 

= 2 



, where we have defined 



U = |dVdV |0(r)p^-^|0(r')P, 



(2.12) 



(2.13) 



pjiO 1^ 



(2.14) 



for q 7^ 0. It is obvious that the on-site Coulomb energy 
U depends on the profile of the Wannier function, which 
expresses the charge distribution near the atomic site. 
Here, we put fq=o = 0, taking into account the charge 
neutrality in the present jellium model. We call the on- 
site U and the inter-site Vq the short-range and long- 
range parts of the Coulomb interaction, respectively. If 
we put Vq = 0, the present model is reduced to the pure 
Hubbard model on a lattice, while if we put ?7 = 0, it is 
reduced to the electron gas model in the continuum space 
except that the momenta are replaced with the crystal 
momenta. 

It is convenient to write eq. (2.4) with matrices as 
Vqd-Q + {Vq + U)ai 



= ^ 

qaa' 

1 



— ^ ( riqi nq^ ) Vq 



Uql 



(2.15) 



where we have introduced the Pauli matrices and the 
unit matrix, cti, (T2, 1T3, and ao, respectively. Here, we 
have defined nqo- = J2k 4ca^i'+qa and 



Va 



•'9 
Vo + U 



u 



We could also rewrite eq. (2.4) as 
1 v-^rl, U, 



(2.16) 



(2.17) 



with 



nq] + Uqi and = (rig^ - nqi)/2. 



2.2 Approximation for the long-range part 

Now, we examine the long-range parts Vq in detail. For 
small q 7^ 0, we could approximate Vq by 



d-y 

VceU VcsuIq 



2 ■ 



(2.18) 



Here, we note that a periodicity relation Vqj^K = Vq holds 
in the original expression eq. (2.14), where K denotes 
any reciprocal lattice vector, whereas it does not hold in 
the approximate expression eq. (2.18). Therefore, if one 
uses eq. (2.18) in calculating the gap function A(fc), the 
reciprocal lattice periodicity A(fe) = A(fe -|- K) of the 
gap function is broken. 

Therefore, we note that the original expression of Wg, 
eq. (2.14), is large near q = K, not only near q = 0, for 
the lattice periodicity. In the gap equation, the momen- 
tum q corresponds to the momentum transfer k — k' , 
where fe and k' denote the electron momenta near the 
Fermi surface. Even when k and k' are in the first 
Brillouin zone, q — K = k — k' — K can be small, 
for example when the Fermi surface is open. For such 
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Fig. 1. Estimation of the nearest neighbor Coulomb interaction 
Va calculated with the approximate equation (2.19) in the finite 
size system with the linear dimension Nsx ■ The broken line is the 
guide for eyes. The thin solid line shows the exact value Va = 1- 



k and k' , we should use the approximate expression 
Vq = Awe'^ /Vceii\q — -K'P, which is very different from 
eq. (2.18) near q = K. Therefore, a more appropriate 
expression of eq. (2.14) is 



max 

K 



47re^ 



V,,ii\q-K\ 



(2.19) 



rather than eq. (2.18). This expression is appropriate not 
only near g = 0, but also near q = K. Near the Brillouin 
zone boundary of q. the derivative of the right-hand-side 
of eq. (2.19) jumps, but since Vq is very small there, it 
does not cause any difficulty in practice. 

Let us examine the error of the above approximation. 
Since eq. (2.19) becomes exact relation for long wave- 
lengths, it is expected that the error of eq. (2.19) is 
larger for shorter wavelengths. Hence, the nearest neigh- 
bor Coulomb interaction 



Va = N ^ ^ Vqe 



(2.20) 



would have the largest error, where a = (a, 0,0) de- 
notes one of the tmit lattice vectors, since the on-site 
Coulomb interaction U is not included in the definition 
of Vq, eq. (2.14). Figure 1 shows the result of eq. (2.20) 
calculated with the approximate equation (2.19) in the 
electron system on the cubic lattice with the number of 
the lattice sites N = N^^, where Ngx denotes the num- 
ber of the lattice sites in one direction. It is found that 
as the system size increases, the value of Va estimated by 
eq. (2.19) approach a value w 1.05 x e^/a, which is only 
5% larger than the exact value. Therefore, the present ap- 
proximation does not cost major error even in the short 
wavelength behavior. 



3. Coulomb Screening 

In this section, we examine the screening effects on 
the Coulomb interactions between two electrons, the 
electron-phonon interactions, and the phonon propaga- 
tors. For convenicuice, we define appropriately scaled 
functions of the charge and spin fiuctuations in the pure 
Hubbard model i.e., the model with U ^ and Vq = 0, 
in the RPA by 



where we have defined 



2X0 (g) 



Uxoiq) 

xo{q) 



1 
1 

21 -Uxoiq)' 



(3.1) 



fc 



G(o)(fc + <z)GW(fc) 



k 



ii^m + £,k+q — Cfe 



(3.2) 



with the unperturbed Green's fiuiction CJ"\ Here, we 
have introduced the four-momentum notation such as 
q = {q,ii'm) and k = (fe,ia;„). 

3.1 Fluctuations and Screened Coulomb Interactions 

In this subsection, we derive an expression of the 
screened Coulomb interactions in the presence of both 
the short and long-range parts of the Coulomb interac- 
tions, U and Vq. 

First, we consider the ring diagrams by the RPA as 
depicted in Fig. 3, where the momenta and spins are 
assigned as shown in Fig. 2. Since the dashed lines of 
the Coulomb interactions have the same momentum q 
in those diagrams, they become divergently large at the 
same time for small |<7|. Therefore, this series gives a 
major contribution for long wavelengths. 

The vertex due to the ring diagrams is obtained as 



r'-s{k,k',q) = Vq ao+Xo{q)% 



(3.3) 



in the matrix form in the spin space. Equation (3.3) is 

rewritten as 



^"'^'^'''''''''^-ll + {2vq + U)X0{q) 



1 1 
1 1 



1 



-U 



1 -1 

21 -Uxoiq) y-'^ 1 
and we could put the elements in the form 

rZf{k,k',q)=TM+Triq)aa' 



(3.4) 



(3.5) 



In the static approximation = 0, we have the effective 
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interaction by the ring diagrams 



q 



U/2 



+ 2K + [7/2)xo(q,0) 



U 



nail- 



(3.6) 



1 



Here, we note that the rotational symmetry in the spin 
space is not retained if we take only the ring diagrams. 

In order to recover the rotational invariancc, wc need 
to take the particle-hole ladder diagrams depicted in 
Fig. 4, with respect to U. The second term in the braket 
of cq. (3.6) becomes divergently large for the momentum 
q near the nesting vector Q, when the Fermi surface 
has a good nesting condition, as known in the Hubbard 
model."''''' It is also known that the contribution from 
the ladder diagrams in Fig. 4 is of the same order as that 
from the ring diagrams. 

In contrast, in the RPA in the pure electron gas model 
with U = 0, the ladder diagrams are ignored. Since the 
momentum of the interaction lines do not always coin- 
cide as those in the ring diagrams, the contribution from 
the ladder diagrams is smaller than that from the ring 
diagrams. Hence, we omit Vq in the ladder diagrams also 
in the present model. 

Thus, the contribution from the ladder diagrams is 
obtained as 

^ (l-(5aa')- (3.7) 



1-Uxoik-k' ~q) 
In the static approximation, wc have 



(3.8) 



where we have defined 



and Sg = {S^)K Adding eqs. (3.6) and (3.8), we obtain 
a rotationally invariant effective interaction 

Vg + U/2 



He 



-T\-- 



+ 2{vq + U/2)x,{q,Q) " " 

-Sq ■ S-q — Unq-^n-ql 



(3.9) 



1-Uxo{q,0) 

where the last term subtracts the double counted con- 
tribution of the first order diagrams in Figs. 3 and 4. 
From eq. (3.9), one could reproduce the vertex parts by 
the RPA both in the electron gas model and in the Hub- 
bard model by taking appropriate limits, U ^ and 
Vq ^ 0, respectively. When we use the effective Hamil- 
tonian of cq. (3.9) in the electron self-energy, we should 
subtract another double counted contribution from the 
second terms in Figs. 3 and 4.^^ 

Within the same approximation, the propagators for 



the charge and spin fluctuations are obtained as 

2X0(9) 



Xciq) 



1 + (2^;, + U)xo{q) 



^^^«)^2i-[/xo(g) 

Hence, we can rewrite eq. (3.9) as 

He = He + Hcf + Hsf, 

where 

^rf = -2^EK + ^)'xc(g,o)r 



'qlV — q 



H.i = -^Y.^U\M,^)Sq-S.q. 



(3.10) 



(3.11) 



(3.12) 



The terms H^i and Hsi can be regarded as the inter- 
actions mediated by the charge and spin fluctuations, 
respectively. 

3.2 Dielectric function and the screening length 

The first terms of eqs. (3.4) and (3.9) are the interac- 
tions between the charge degrees of freedom of electrons. 
Hence, we obtain an expression of the dielectric function 



K{q) = l + 2{vq + %)xo{q) 



(3.13) 



In the long wavelength limit q — > 0, since Vq ^ U, we 
obtain Thomas-Fermi dielectric constant 



K{q) = 1 + ^ = KTF(g) 

q 



(3.14) 



with gxF^ = 87re^A^(0)/KoU, where Atf = 1/9tf and 
A^(0) denote the; Thomas-Fermi screening length and the 
electron density of states at the Fermi energy per a spin 
and a site. For the short wavelength \q\ » qg, we can 
omit Vq compared to U, 



K{q) = 1 + Uxo{q) = Ku{q), 

which is the result in the Hubbard model. 
Equation (3.13) is written in the form, 

K{q) = K,u{q)Kv{q), 

where 



(3.15) 



(3.16) 



K,{q) = l + VqXi°\q) = l + 



91 



Qa — QtF 



2N{0) Ku{q) 



(3.17) 



2_ 2 X0(?) 

Is — Itf 



N{0) 

with Xc^^ defined in eq. (3.1). In the absence of U, in 
the system with the spherical Fermi surface, eq. (3.16) 
is reduced to the Lindhard's result of the RPA in the 
electron gas model. 
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k-q,<J 
\ 



/ 



k'+q,a' 
/ 



\ 
k',a' 



Fig. 2. Assignment of the momenta and spins in the two body 
vertex part r(fe, k' ,q). 



Fig. 4. Summation of tlic particlc-liolc ladder diagrams. The def- 
initions of the lines are the same as those in Fig. 3. The incom- 
ing and outgoing thin solid lines are to show how to apply the 
diagrams to the two-particle vertex in the self-energy and the 
superconductive gap equation. 



Fig. 3. Summation of the ring diagrams. The solid and broken 
lines denote the electron Green's function and the Coulomb in- 
teraction eq. (2.12). 





For long wavelengths and in the static approximation, 

we have xo{q) » N{0) and Ku{q) w 1 UN{0). Hence, 
we obtain the screening length 



(3.18) 



Ascr = Atf ^/ku{0) = Atf VT+Jtc 
from eq. (3.17), where we have put 

/ic = UN{0). (3.19) 

Therefore, it is found that the screening length is length- 
ened by a correlation effect due to the on-site Coulomb 
repulsion in comparison to the Thomas-Fermi screening 
length Atf ■ It is plausible that the local electron density 
less easily deviates from the uniform value in the presence 
of strong correlations. In other words, the charge fluctu- 
ations are reduced {[xi^^]u^Q < [xi^hu=o) and ku > 1, 
when > 0. It is interesting that the short-range part of 
the Coulomb interactions affects the long-range behavior 
through a many body effect. 

3.3 Corrections to the electron- phonon interactions 

Now, let us examine renormalization of the electron- 
phonon interaction due to the Coulomb screening effect 
as depicted in Fig. 5. As the screened Coulomb interac- 
tion, we take the series of the ring diagrams shown in 
Fig. 3. However, if we apply the ladder diagrams in thick 
dashed line in Fig. 5, the ring part consists of many elec- 
tron lines. In such a ring, the integrand could be large 
only in a very small part of the phase space of the in- 



Fig. 5. Renormalization of the electron-phonon interaction. The 
solid lines denote the electron Green's functions, and the closed 
circle denotes the bare electron-phonon interaction. The thin 
solid and wavy lines are to show how to apply the diagrams 
to the electron and phonon Green's functions, respectively. The 
thick dashed line denotes the screened Coulomb intera<;tion. 



tegral variables. Hence, we omit the ladder diagrams for 
simplicity. This simplification is also justified from the 
rotational invariance in the spin space. Since we could 
take the spin summation for the ring part in Fig. 5, the 
spin interaction in eq. (3.4) vanishes. Thus, the electron- 
phonon vertex part is not rotationally invariant if we 
retain the ladder diagrams. Therefore, we ignore them, 
and obtain the electron-phonon vertex part 



M{q) 



(3.20) 



3.4 Corrections to the phonons 

In this subsection, we examine the phonon Green's 
function. We consider the renormalization of the phonons 
due to the Coulomb screening through the electron- 
phonon interaction shown in Fig. 5. Then, we obtain the 
diagram equation shown in Fig. 6. This approximation 
is equivalent to that for the polarization propagator as 
shown in Fig. 7. 
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VWV = ^AAA/ + 



Fig. 6. The diagram equation for tlie dressed phonon Green's 
function. The thiclc and thin wavy lines denote the dressed and 
bare phonon Green's functions, respectively. The renormalized 
electron-phonon interaction is shown in Fig. 5. 



Fig. 7. The diagram equation for the polarization propagator. 
The solid line and thick dashed line denote the bare electron 
Green's function and the screened Coulomb interaction. 



The phonon Green's function is defined by 

V{q,T) = -{Tr[ip^iT)4]) (3.21) 

with (fq = bq + whoro we have defined A{t) = 

e^^ Ae~'^^ . The bare phonon Green's function 2?o(g, r) 
is defined in a similar manner, and the Fourier transform 
is expressed as 

Vo{q,iu^) = -^-fr^. (3.22) 

The polarization propagator of Fig. 7 summed over 
spins is expressed as 



Piq) = 2xo(9) - [xo{q)? J2(^c{q) + Tf (<z)c7a') 



(3.23) 



2xo(g)/K(<?) = xdq), 



where Xc{q) denotes the charge fluctuations. From Fig. 6 
with Fig. 5 and eq. (3.20) or with Fig. 7 and eq. (3.23), 
we obtain 

2M?Yo(g) 

•Diq) = Vo{q) - Po(g) 7, ^ P(g), (3.24) 
and thus 

V-'=V^'+M'^Xc{q). 
Using eqs. (3.22) and (2.5), we obtain 



(3.25) 



■D-' = -^K + Mq)f 



with 



u){q) = 



' Ku{q) _ 



(3.26) 



(3.27) 



i^iq) iiviq) 

The renormalized phonon dispersion energy is obtained 



by solving Wq = oj{q,Cjq), where the function ui{q,ui) is 
obtained by the analitic continuation co + iS from 

oj{q,'Wm)- 

In the long wavelength limit q — > 0, we obtain 



U)q = c\q\ 



(3.28) 



with c = flp/qg. As is widely known, the longitudinal 
optical phonons are screened by electrons into the ac- 
coustic mode with a dispersion energy proportional to 
\q\ as expressed in eq. (3.28). It is easily verified that 
putting U = 0, one could recover all the standard results 
of the RPA of the electron gas model in the back-ground 
jellium positive charge.^ For U ^ 0, we find that the 
electron correlations due to the on-site Coulomb repul- 
sion weaken the screening effect and increase the sound 
velocity. 

4. Superconductivity 

In this section, we examine the two-particle vertex part 
r„cr^{k,k'), which induces superconductivity. We clar- 
ify the corrections to the vertex part due to the on-site 
Coulomb interactions. 

4-1 Definition of the effective interactions 

The two-particle vertex part Fo-cn {k, k') is defined as 
one which scatters the pair of two electrons with {k, a) 
and {—k,ai) to that with {k',a) and {—k',ai). The gap 
equation of superconductivity is written as 

(4-1) 

X G,{k')C{-k')A,,,{k'), 



where we have introduced an abbreviation 



(4.2) 



and Ga{k) denotes the dressed electron Green's function. 
The vertex part T„„-^ plays a role of the effective inter- 
actions. The gap equation (4.1) includes the dynamical 
effect by the Matsubara frequency. 

The electron self-energy is calculated by 

S4fc) = -^r('^)(A;,fc')G.'(fc'), (4.3) 

where F^'^] denotes the appropriate vertex part for the 
self-energy consistent with Tuu-^ . The consistency is con- 
sidered phenomenologically or microscopically with the 
Ward-Takahashi identity. Equations (4.1) and (4.3) are 
called Eliashberg equations. 

In order to estimate accurately, we need to solve 
the full set of equations. Especially when we consider 
the spin-fluctuation-mediated superconductivity, since 
the attractive part of the interactions is much weaker 
than the repulsive part, the self-energy effect is impor- 
tant.'^^ However, since in this paper we are mainly in- 
terested in the phonon-mediated pairing interactions, we 
omit the self-energy renormalization of the electrons from 
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now on. 

We can improve the perturbation theory^^'^^ by in- 
cluding rcnormalization of the electron Green's function 
and the vertex corrections in Ta^criikjk'). For example, 
we could take into account the correction to the spin 
and charge fluctuations from the self-energy effect by re- 
placdng xo{q) with xil)^ as in the renormalized RPA, 
although it requires heavy numerical calculations. 

However, we expect that the most essential diagrams 
are taken into account in the RPA. Within the RPA, 
we could reproduce a physical situation that the strong 
charge and spin fluctuations mediate the pairing inter- 
actions. This consideration is partially based on a phe- 
nomelogical viewpoint. However, more or less, we could 
not avoid such a phenomenological consideration in the 
perturbation theories based on effective Hamiltonians, 
such as the Hubbard model. In any perturbation theory, 
we need to ignore almost all diagrams in the infinite se- 
ries of the perturbation expansion. Taking a particular 
series of the diagrams can be rigorously justified only in 
an appropriate limit, but otherwise it is not quantita- 
tively accurate very much. Furthermore, we need to use 
the effective Hamiltonian, which is usually far simpler 
than the real materials. Usually, effective Hamiltonians 
are useful for qualitative study, but are not very precise 
quantitatively. Therefore, even if we develope a pertur- 
bation theory by purely theoretical efforts beyond the 
RPA, the applicability of the quantitative results to the 
real materials are limited, unless any qualitatively dis- 
tinct result is obtained. 

4-2 Expression of the effective pairing interactions 

In this subsection, we derive the expression of 
^aai{k,k') from the results of the previous section §3. 
From the momentum and spin assignment depicted in 
Fig. 2 that r™?(fc,-fc,fc - k') and F^^^^ ^ _ ^/^ 

in eqs. (3.5) and (3.7) contribute to the vertex part 
PcrcTi {k, k'). 

The contribution from the charge fluctuations and the 
screened Coulomb interaction is obtained from eqs. (3.4) 
and (3.5) as 



rc(9) = T.: 



2Va + U 



U/2 



21 + {2vg + U)xo{q) Kiq) ' ^^'^^ 

where we have put q = k — k'. Similarly, the contribu- 
tion from the spin fluctuations with respect to the Sz 
component is written as 

U 



(4.5) 



2 l-C/xo(g)^ 

In addition to them, we obtain the effective interaction 
mediated by phonons 



rph(9) = mq)fV{q) = 



q 



-2Vtp/n{q) 



(4.6) 



K{q) i^l + [Cj{q)? 

with Cj{q) given in eq. (3.27). 

The contribution from the spin fluctuations with re- 



spect to the Sx and Sy components is obtained from the 
ladder diagram as 

^^crcriPg" ) = (5(j(j^r (A;, A;, k ) 

^ U (4.7) 

^'^'^^l-Ux^ik + k'Y 

where 5aa^ = 1 — 6aai ■ By the parity of the gap function 
A(fc,ia;„), we could replace the vertex function eq. (4.7) 
in the gap equation eq. (4.1) with 



f^CTCTiP^ {k,k ) — sSaai 



u 



where we take the sign s = +1 and s = — 1 for singlet 
pairing and triplet pairing, respectively. 

Collecting the contributions of eqs. (4.4) - (4.8), we 
obtain the vertex part 

T^,, {k, k') = r,{k - k') + aaiTfik - k') 

+ T^^{k-k') + ~5^^,Vt-{k,k') (4.9) 

- U Jo-o-i . 

Here, the last term —Ud^fr^ is to subtract the double 
counted contribution as mentioned below eq. (3.9). With 
a definition 

eq. (4.9) can be rewriten as 



= Fc - CTCTiFs + Fph + 25„„^sT^ -US. 



(4.11) 



Therefore, for singlet pairing, the vertex part Tcr,ai {k, k') 
is expressed as 



■ SP 



Fc + 3F, + F 



ph • 



while for triplet pairing 



-tp 



Tc + Fs -I- Fph. 



(4.12) 



(4.13) 



If we put Vq = 0, the spin and charge fluctuation terms 
are reduced to the expressions obtained in many arti- 
cles. Due to the complicated momentum de- 
pendence, it can work as an attractive interaction or as 
a repulsive interaction in the formation of the Cooper 
pairs, depending on the symmetry of the gap function 
and the shape of the Fermi surface. 

Since we have obtained the expressions of the vertex 
part eqs. (4.12) and (4.13), it is possible in principle to 
solve the Eliashberg equations (4.1) and (4.3). However, 
we leave it for future study, and derive weak coupling 
effective Hamiltonian in the next section §5. 

4-3 Physical interpretation 

In this subsection, we rewrite the expressions obtained 
above for a physical interpretation. We also discuss the 
expressions in the two limits, C/ ^ and Vq 0. 

In eq. (4.11), the screened Coulomb interactions and 
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the phonon-mediated interactions are rewritten as 



ph 



Vg + U/2 



u 

2ku 



KKu -w^ + [oj{qW 

,2 



(4.14) 



UJ 



KKu — [li'(g)]^ ' 

where we have made the analytic continuation ivm 
u! ± i5, and put q = k — k'. 

Therefore, the total two-particle vertex part is given 
by 

(4.15) 



r<,.,(g) = r[,°) + 



KKu — [w(q')]^ 

with q = k—k', where the function r^*^ in the first term is 

nothing but the vertex part in the pure Hubbard model, 
i.e., the model with [/ ^ and Vg = 0.22.35,36 p^j. gingiet 
pairing, it is expressed as 

1 U 3 U 



r(o) - 
'-u — 



2 1 + (7x0 
U 

~ 1 - U^x'o 
which can be written as 



+ 



2 1 - C/Xo 

U 



U 



(4.16) 



U, 



u 



lu\(^^+3U\r^ + U (4.17) 



with eq. (3.1). For triplet pairing, it is expressed as 



1 



U 



+ 



1 



u 



r(o) ^ 

^ 2 1 + C/xo ' 21-C/xo 



U 



u^xV 



which can be written as 



(4.18) 



(4.19) 



In eqs. (4.17) and (4.19), the physical interpretation of 
each term is clear.22>35,36 r^j^g gj.g^ terms of eqs. (4.17) 
and (4.19) arc the interactions mediated by the charge 
fluctuations, while the second terms are those by the 
spin fluctuations. Obviously, the last constant terms in 
eqs. (4.17) and (4.19) express the bare on-site Coulomb 
repulsive interactions, and do not contribute to the gap 
equation for anisotropic pairing. If we put Vg = 0, we 
obtain 



.(0) 
■ U ' 



(4.20) 



and all the equations are reduced to those in the Hubbard 
model. 

The second term of eq. (4.15) corresponds to the sum- 
mation of the screened Coulomb and phonon-mediated 
interactions. If we put J7 = 0, we obtain 

,2 



CO 



K„ U)^ 



(4.21) 



where Ky = 1 + 2vgXo and uj{q) = ftp/ ^yKy{q). This 



equation is easily rewritten as 

T,aAQ) = - + [M{q)?V{q) 



Ky 



V. 



M2 

q_ _l_ q 



(4.22) 



Ky U!^ - [w(g)]2' 

We note that the expression of T^^ri in the presence of 
both U and fg, i.e., eq. (4.15), is not obtained by sim- 
ply adding eqs. (4.20) and (4.21). The second term of 
eq. (4.15) includes the corrections due to C/ in a rather 
complicated manner as explained in the next subsection. 

4-4 Interpretation of the corrections due to U 

Now, let us interpret how the corrections from U en- 
ters in the phonon-mediated pairing interaction, i.e., the 
second term of eq. (4.15). We can derive eq. (4.15) in 
another manner as follows, starting from eq. (4.22), i.e., 
the effective interaction in the absence of U. 
First, we rewrite eq. (4.22) as 



2Ml 



Ky KyVtp IjP' — [W(g)]2 ' 



(4.23) 



and in this equation wc replace Ky with ky rcnormal- 
ized by the Hubbard-type charge fluctuations as shown 
in eq. (3.17). At the same time, we replace the phonon 
frequency Lij{q) with cj(g) of the fully dressed phonons, 
which is given by eq. (3.27). Secondly, we modify the 
electron-phonon coupling constant Mq with 

^ = M,[l-[/xo + C/^xgT---]. 

Ku 

The series in the right-hand-side is easily verified from 
the diagraniatical technique. At the same time, we insert 
\Iku into the both end of the screened Coulomb line 
which corresponds to Vg/iiy^q). In other words, we re- 
place Vg/ky with Vg/kyKij. Thirdly, we add the Hubbard- 
type fluctuation diagrams T^^ including bare U. Then, 
we obtain the expression of T^ai (?) 



r^^i (q) = T 



(0) 



2M2 



2 ~ 



which coincides with eq. (4.15), since 



1 v„ 



2Ml 



(4.24) 



(4.25) 



holds from eqs. (3.16) and (2.5). 

5. Effective Interactions for Weak Coupling 

In this section, we derive the effective pairing interac- 
tions along the scheme of the traditional weak coupling 
theory. The effective interactions derived here are not 
suitable for accurate estimations for the transition tem- 
perature, but useful for qualitative and semi-quantitative 
arguments. For quantitative purpose, direct numerical 
calculations based on the effective vertex eqs. (4.12) and 
(4.13) is more suitable. 
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5.1 The case of U = 

First, we briefly review the standard derivation of the 
efFective pairing interactions for U = 0.^ In this case, the 
total two-particle vertex part Tg-a--^ is given by eq. (4.21). 
In this equation, it is explicitly shown that the phonon- 
mediated interaction is attractive for small lu overcoming 
repulsive Coulomb interactions. This effect is called the 
overscrecning effect. 

In the weak coupling theory, the frequency dependence 
of the two-particle vertex part T(q) is simplified by a 
replacement with a step function as shown in Fig. 8, that 



is, 



r{q,u;) 



[l-a0{u,-\oj\)]. 



(5.1) 



K{q,0) 

where ujc and a denote effective constants which express 
the frequency range and the depth of the attractive part 
of the interactions. Here, the static approximation for the 
dielectric function has been introduced. 

The effective constants a and ojc are complicated quan- 
tities that reflect the dispersion energy and the density of 
states of the dressed phonons. However, it is reasonable 
to consider that the constant uic is of the order of the 
Debye frequency wd, i-e., the upper limit of the dressed 
(observed) phonon frequency cuq. According to the stan- 
dard weak coupling theory, we regard those constants as 
parameters in this paper. 

For isotropic superconductors, the prefactor Vq/ K{q, 0) 
is often replaced by an effective constant Vq. We should 
note that Vb is an averaged quantity of Vq and different 
from the on-site Coulomb interaction U derived in §2.^' 

Further, in the weak coupling theory, the cutoff lOc 
in oj integral is taken into account by introducing the 
same cutoff in ^ integral, where ^ denotes the electron 
energy measured from the Fermi energy. In other words, 
we replace the cutoff function 9{(jJc — I'^l) with d{oJc — 

laiM^c-l^fc'l). 

After these simplifications, wc could solve the gap 
equation along the scheme of the standard weak cou- 
pling theory. The Coulomb repulsion, i.e., the first term 
in the braket of eq. (5.1), is taken into account through 
the effective Coulomb parameter defined by 



(5.2) 



1 + VoN{0) In(I^cK) ' 

where Wc denotes the cutoff energy of the Coulomb 
interactions, which is of the order of the band width 
W .^^ The influence of the Coulomb repulsion is weaken 
(Mvo < K)-^(O)) by the retardation of the phonon prop- 
agation when LOc < Wc- The superconducting transition 
temperature Tc is given by 

Tc = l.lSwce"^/^^"''^''^ (5.3) 

where A = aVoN{0). 




Fig. 8. Schematic figure of the two-particle vertex part r(g) as a 
function of ui/ojc. The solid, dashed, dot-dashed lines show r(g) 
as functions of u> with different fixed values of q. The thick dotted 
line shows a step function which is introduced to approximate 
r{q) in eq. (5.1). 



5.2 The case of U ^ Q 

We extend the above formulation to the case of [/ 7^ 0. 
We introduce the same simplification as eq. (5.1) in the 
context of the weak coupling theory. Hence, it is obtained 

that 



r„,,(q,c^) = r[,"'+i)q [1 



(5.4) 



with Vq = Vq/ K{q,0)Kij{q,0), where the static approxi- 
mation has been also introduced in k, kij and r(,°) as m 
eq. (5.1). 

Let us compare the expressions for U = Q and f/ 7^ 0, 
i.e., eqs. (5.1) and (5.4). One of the important differences 
is that eq. (5.4) includes the Hubbard-type effective in- 
teraction r[j'' that wc have discussed in §4.3. 

Another difference between eqs. (5.1) and (5.4) is the 
renormalization of the second term of eq. (5.4) due to the 
on-site Coulomb interaction U. In this term, the prefac- 
tor Vq and the constants a and lUc are modified. The con- 
stants a and cOc are usually treated as parameters in the 
weak coupling theory, and sometimes determined by ex- 
periments. Hence, we do not need to know the differences 
from their original values, and the modifications in a and 
LOc do not need to be considered explicitly. Therefore, the 
important modification is in the prefactor Vq, because it 
affects the momentum dependence. Prom eqs. (3.16) and 
(3.17), it is written as 



Vq = ^- 

KKJJ 



1 



91 



Putting U = in. eq. (5.5), it reduces to the factor of 
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eq. (5.1), i.e.. 



Hence, we write it in the form 



2X0 \q\- 



It 



(5.6) 



Comparing eqs. (5.5) and (5.6), we note that is re- 
placed with a smaUer quantity qg in eq. (5.5). As dis- 
cussed in §3.2, the screening length Agcr = 9tf = [9s]g=o 



is lengthened into Ag 



because f/ > 0. 



Then, the peak of eq. (5.5) becomes sharper, and the 
anisotropic components increase. 

In summary of this part, we find two effects of the 



(0) 
u ' 



short-range correlations: (1) The Hubbard term T 
which includes bare U, is added as pointed out in the pre- 
vious section §4; (2) The short-range correlations change 
the momentum dependence of the phonon-mediated pair- 
ing interaction. Both effects unfavor isotropic pairing, 
but favor anisotropic pairing of some symmetries. 

It is also found that for this mechanism of triplet pair- 
ing, the existence of the on-site Coulomb interaction U is 
essential, but not the s-wavc (isotropic) component of Vq, 
such as Vq. We note in the second term of eq. (5.4) that 
the long-range part Vq of the Coulomb interaction and 
the phonon-mcdiatcd pairing interaction —avq9{ujc— 
have the same momentum dependence. Hence, it is not 
effective for the change in the dominancy of the coupling 
constants. In fact, when wc put the dominant isotropic 
component Vq and the subdominant anisotropic compo- 
nent V\ , the difference of the effective coupling constants 
which include the negative contributions of the effective 
Coulomb parameters is calculated as 



1 + Vqx 
= {Vo - Vi) 



aVi 



1 + Vix 



1 



a — 



(i + yo3:)(i + yi.T) 



(5.7) 



where x = N {0) \n{Wc / u^c) ■ We have a > 1 when the 
overscreening effect gives rise to a net attractive interac- 
tion for low frequencies. Thus, when Vq > Vi, the right- 
hand-side of eq. (5.7) is positive. Therefore, it is found 
that the long-range part Vq, which gives rise to the third 
and fourth terms in the left-hand-side of eq. (5.7), does 
not change the sign of eq. (5.7). In the rest of this pa- 
per, we investigate the change in the dominancy of the 
coupling constants due to the Coulomb interaction more 
explicitly. 

5.3 Interactions m,ediated by spin fluctuations 

In this subsection, we develope a weak coupling the- 
ory including interactions mediated by spin fluctuations 
through r[^^ In particular, when the spin fluctuations 

are relevant, the vertex part T^\q,u>) consists of two 
parts as a function of q and u: (a) A broad peak that 
reflects the energy scale of the band width ~ W; (b) 
A sharp peak that reflects the energy scale of the spin 
fluctuations cogi, as discussed in our previous papers. ^^'^^ 



rPiq,uj) « f'^>{q,0) + T':^>{q,0)9{co,{ - \u;\), (5.8) 

where the first and second terms which correspond to 
(a) and (b) mentioned above, respectively. Obviously, 
r[^'' includes the on-site Coulomb repulsive interaction. 
Therefore, we obtain 

T^^,iq,uj) = Vc{q) + T'^j^\q,0)eiujsf 
- avq 0{loc -\oj\) 

from eq. (5.4), where we have defined a Coulomb term 

Vc{q) = t'uHQ,0) + ig. 

As mentioned in §5.1, we replace 0{LVk — with 
d{u)k — \^k\)0{uJk — l^fc'l) within the weak coupling the- 
ory, where Ui = lOc and ui2 = i^ai- Therefore, we obtain 
effective interactions 

Kff (fc, k') = Vcik- k') + V,{{k, k') + yph(fe, fc'), (5.10) 

where 

y,f(fc,fc') =r(°)(fc-fe',0)^(u;,f- |ai)e(a;sf- ICfe'l) 
Vph(fc,fe') = -avk-k'0{uic - \^k\)9{oJc - l^fe'D- 

(5.11) 

When u>s{ <C W and Wc <C M^, we may regard F^f ' and 
Vk-k' as functions of fc|| and fcy, which denote the mo- 
mentum components of k and k' along the Fermi surface. 
Hence, we introduce a complete set of the basis func- 
tions 7a(^||) with a = 0,1,2,---, which are normalized 
by the average on the Fermi surface. For example, we 
may replace fcy with k = k/\k\ or {9,(p), and take a = 
s,Px,Py,Px ■ ■ ■ or a = {l,m) in the spherically symmetric 
systems. We define V;f(fc||,fc||) = [r^f^(A; - fe', 0)]|^=^j|^=o 
and yph(fc||, fc||) = — Q:[{ife_fe']|^=|j^=0) and expand the in- 
teractions as 



=,(0) 



(0)/ 



(5.9) 



Vfe(/C||,/Cf|)=^yfc"7a(fc||)7a(fc||), 



(5.12) 



where the index k donotes the kind of interactions, for 
example, Vi = Vph, V2 = V^f and V3 = Vc- We define 
dimensionless coupling constants Apj^ = — V^^A^(O) and 
A"f = —V^N{0), and cflFcctive Coulomb parameter /x* = 
VSNiO)/[l + VSN{0) In(T^cM)]. 
The gap function is also expanded as 



A(fe) = ^A„7„(fc||) 



(5.13) 



on the Fermi surface. By taking the basis functions 
appropriately, we can make the gap equation diagonal 
with respect to the symmetry index a upto the linear or- 
der of A. While each decoupled gap equation gives a dif- 
ferent second order transition temperature T^a, we must 
take the highest one of Tea's, since the system exhibits 
superconductivity when at least one of A^'s is nonzero. 

In the same way as in our previous paper, we can ob- 
tain the transition temperature Tc and the isotope eSect 



J. Phys. Soc. Jpn. 



Full Paper 



Hiroshi Shimahara 13 



coefficient am defined by 



^InTc 



where M denotes the relevant atomic mass. Because the 
application is straightforward, we shall omit the deriva- 
tion and only present the results. 

When LUgi < Wc ~ wd, which occurs in the prox- 
imity of the magnetic instability, we obtain the transi- 
tion temperature Tc = 1.13cJcexp[-l/(A^f + A^^)] with 



A: 



ph 



(A"h - - (A"h - iil)Hco,/io,f)] and the 



isotope effect coefficient 



"IE 



VA«f + A«jj; L VA«,-M*>' 



(5.15) 



When LOsi > llJc ^ ^^D , which occurs away from the mag- 
netic instability, we obtain the transition temperature 
Te = 1.13a;cexp[-l/(A«h + A"/)] with A«; = (A«f - 
MS*)/[l-(As"r-MsT)lnKfM)] and m,T = ^c"^(0)/[l + 
VqN{0) ln{Wc/tUsi)]- The isotope effect coefficient is ob- 
tained as 



1 



aiE 



^ ) 



2i 



(5.16) 



5.4 Anisotropic pairing mediated by phonons 

In this subsection, wc derive expressions of the effective 
coupling constants and the effective Coulomb parameters 
for anisotropic pairing. We concentrate our attention on 
the interactions mediated by phonons, and ignore the 
attractive interactions by F^''. 

Hence, we consider the case in which we could ignore 
the momentum and frequency dependences of Xo('?), so 
that xo{<l) « Xo(0,0) = iV(oj holds. Although this sim- 
plification becomes qualitatively better away from the 
magnetic instability, even there the spin fluctuations still 
have an important effect of enhancement of the repulsive 
Coulomb interactions as 

j-ii-^ = c; + c;',<">(o)>t/. (5.17) 

Further, since T^^^ becomes constant in the simplifi- 
cation mentioned above, it affects only s-wave pairing. 
Hence, we use the expression of Fj^' for singlet pairing 
eq. (4.16). Prom eqs. (5.4) and (5.5), we obtain 



F,,,(q,w)7V(0) 



1-Mc 1-Mc 2(1 -h Ate) 



1 + 



(5.18) 



where we have defined a parameter g hy a = 2gN{Q) 
and 



(5.19) 



with 



91 



9s 



9tf 



(5.20) 



We expand F{q) and A(A;) on the Fermi surface with re- 
spect to the momentum components parallel to the Fermi 
surface as in eqs. (5.12) and (5.13). We define the expan- 
sion factors Ca of F by 

F(fc-fe') =EC'a7a(fc||)7a(fc||)- (5-21) 
a 

For convenience, we write a = for s-wave pairing, while 

a 7^ for anisotropic pairing, from now on. 

Therefore, from eq. (5.18) and the replacement of 
6{ijJc — |w|) with 9{uJc — |^fc|)^('^c — ICfe'l) that is adopted 
in §5.1 and §5.3, the a component of the effective inter- 
actions is written in the form 

K.7V(0) = -X^e{u;, - M)0{oj, - iCfe'l) (5.22) 

where Aq, and denote the dimensionless coupling con- 
stants of the phonon-mediated pairing interaction and 

the repulsive Coulomb interaction, for "a-wavc" super- 
conductivity, which are expressed as follows. The con- 
stant Aa is expressed as 

gN{Q) 



Aq. 



1 + /iC 

for any value of a. In contrast, the constant fia has dif- 
ferent forms depending on a 
pairing (a = 0), we obtain 



(5.23) 
IS dif- 

or a ^ 0. For s-wave 



Co 



(5.24) 



i-Mc 2(1 + ^c)' 

while for anisotropic "a-wave" pairing [a ^ 0) 

For the retardation effect, the effective Coulomb param- 
eter /i* becomes smaller than as 

The derivation of eq. (5.26) is explained, for example, in 
Refs. 1 and 17. The superconducting transition temper- 
ature for "a-wave" pairing is estimated from the above 
parameters by the equation 



Tea — 1-13 lJcC 



.-1/Ao 



(5.27) 



where A^ = Xa — IJ-a- As argued below eq. (5.13), the real 
transition temperature is the highest one among all 

T 's 

We obtain the isotope effect coefficient 



1 



(5.28) 



putting A"f = in eq. (5.15). The reverse isotope effect 
aiE < occurs when /U* > A^, i.e., 2 > Aa//i*. 
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From eqs. (5.23) and (5.25), we have A^/Ma = 2gN{0) 
for anisotropic pairing. Hence, for aiE < to occur in 
anisotropic superconductors, the inequality 1 > gN{0) 
needs to hold, because 2 > A^/^* > Xa/lJ^a > 2gN{0). 
However, for such weak coupling, the transition temper- 
ature becomes extremely low so that it is not practically 
observable. If we take into account repulsive contribution 
to fia from r|y\ which modifies eq. (5.25), the reverse 
isotope effect aiE < is possible even in the present 
situation. 

In eq. (5.26), the Coulomb cutoff energy Wc is usually 
of the order of the band width W. More precisely, it is 
defined by 



ln(W 



ep) + Inep 



(5.29) 



as is easily verified from the derivation of eq. (5.26), 
where W and ep denote the band width and the chem- 
ical potential (Fermi energy) measured from the bot- 
tom of the electron band. Equation (5.29) is written 
as Wc = ^/e^^(W'^-eF)■ Hence, it is easily proved that 
Wc < W/2, where the equal sign holds when gp = W/2. 
In particular, when the charge carrier density is very 
small as in semi-metals, that is, ep <C W, we obtain 
Wc « ^/W^ «; W. 

5.5 Dimensionality and anisotropic pairing 

As we could see in eqs. (5.18) - (5.21), the anisotropic 
pairing interactions can dominate when qs is small and 
C„^o is near Cq. Since q^ is proportional to gxp as seen 
in eq. (5.20), small gxp favors anisotropic pairing. In sys- 
tems with isotropic Fermi surfaces, we have 



Hi) 



ej2ki 



(5.30) 



l-cose + ql/2kl as -cose' 

where as = 1 + q^/2kp = 1 + q^pf^^Fi^ + A*c), and 
denotes the angle between k and k' , since q = fc — fe'. 
Thus, it is found that the scaled parameters q^p/2kp 
and fic are the essential parameters in the sense that 
they determine the property of the function F{q). 

First, we discuss the three dimensional system. In sys- 
tems with spherically symmetric Fermi surfaces,^ we ob- 
tain 



2fcp 



1 nih 



7r(37r2) 1/3^1/3 m an nV3 m 2.6[A] 



(5.31) 



where m, rrih, ob and n denote the bare electron mass, 
the band effective mass, the Borh radius Ob = fi? /me^ « 
O.5293A, and the electron (or charge carrier of ±e) den- 
sity per a site, respectively. Therefore, as the lattice con- 
stant a increases, the ratio q^p/2kp increases, which im- 
plies stronger screening. 

We obtain a qualitatively different result in the quasi- 
two-dimensional system. In systems with cylindrically 
symmetric Fermi surfaces, we obtain 



2kp 



1 rrib 



1 mb 

irn m and n m 1.7[A] x d' 



(5.32) 



where a and d denote the lattice constant in the layer 
and the layer interval. As pointed out in our previous pa- 
per, as the layer interval d increases, the ratio q^p/2kp 
decreases, which implies weaker screening, in contrast to 
the result in the three-dimensional systems mentioned 
above. 

6. Examples of the Application 

Purpose of this section is to show examples of the 
application of the effective Hamiltonian derived in the 
previous section. For simplicity, we consider spherical 
symmetric systems. We have already examined a spher- 
ical symmetric system in our previous paper, ^ but the 
present treatment is based on a more basic model. Con- 
sequently, the corrections due to U are taken into account 
in more detail. The transition temperatures obtained in 
this section should be taken as crude estimations, be- 
cause the model is simplified. Nevertheless, the qualita- 
tive results obtained here reflect essential aspects of the 
present mechanism. 

From the symmetry, we take the basis function as 
7a(fe) = Yim{e,v) with a = {l,m). Therefore, 

F{k -k') = J2 Ci [YimiO, </')]* Yimie', /) 

l,m 

= ^(2/ + l)C;PKcOS0fefe') (6.1) 

I 

A{k)=^AimYimie,ip), 
Im 

where 9^^' denotes the angle between fe and k' . Since 
Poiw) = 1 and Pi{w) = w, we obtain from eq. (5.30) 

dH, Q!s — 1 CKs — 1 



Co = 
Ci = 



1 _ as - 
An as — cos 9 2 

dfl {as — 1) cos^ 



In 



1 



1 



4n a^ 

(as - 1) 



cos^ 



(6.2) 





as , 


as 


- 1 




n 


m 










2 


as 


+ 1 





where dQ = sm9d9d(p. Considering the signs of the 
three quantity, Aq, Ai, and Aq — Ai, we obtain the phase 
diagrams as shown in the following subsections. 

6.1 Typical parameter values for ordinary metals 

Figure 9 shows the phase diagram for the set of pa- 
rameter values Wc/uJc = 100 and q^p/2kp = 0.5. The 
formar ratio Wc/i^c = 100 is typical of most metallic 
systems since Wc ~ W and Uc ~ tou- The latter ra- 
tio q^p/2kp = 0.5 implies the Thomas-Fermi screening 
length Atf = l/fcp ^ o/tt, which seems typical of metals, 
too. Because at nc = UcN{0) « 1, the magnetic insta- 
bility Xs ^ 00 should occur, the results are plotted only 
for UN{0) < 1. It is found that the triplet state occurs 
for large U, which confirms one of our previous results.^ 
It is also found that for large g, the s-wave state is more 
favored than the p-wave state. 
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Fig. 9. The phase diagram on the plane of the on-site Coulomb 
energy U and the effective coupling constant g of the phonon- 
mediated pairing interaction, in the spherical symmetric system. 
The parameters are chosen as g^p/2fcp = 0.5 and WqI'jJc = 100. 
The abbreviation SC stands for superconductivity. 
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Fig. 11. The transition temperatures as functions of the effective 
constant g of the phonon-mediated interaction for the narrow 
band case, in which we set q^p/2kp = 0.5, Wq/oJc = 5. We have 
put U!c = 300 K as an example. The solid and dashed lines show 
the transition temperatures Tcp and T^s of the p-wave and s-wave 
states for UN{0) = 0.5. For each value of gN{0), the higher one 
of Tcp and Tcs is the final result of Tc, i.e., the physical Tc. The 
thin solid and dashed lines show those for UN{0) = 0.7. 
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Fig. 10. The phase diagram for a narrow band case, in which we 
set q^p/2kp = 0.5 and Wc/u)c = 5. Along the vertical thin 
dotted line, the transition temperatures are plotted in Fig. 11. 



Therefore, for the present set of parameter vahies, the 
resultant dimensionless coupling constant Xp for the p- 
wave state is quite small, where it overcomes the s-wave 
state. For example, when /ic = UN{0) = 0.5, the high- 
est transition temperature to the p-wave state occurs 
near the phase boundary at gN{0) « 1.05. However, 
at fxc = 0.5 and gN{0) = 1.0, we obtain Ao fa 0.031 
and Ai w 0.037, and such a small value of Ai gives 
Tcp ^ 10^^°K, which is too low to be observed in prac- 
tice. For gN{0) > 1.05, we obtain larger T^., but there 
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Fig. 12. The dimensionless coupling constants Aq, as functions 
of Wc/^c- The solid and dashed lines denote Ap and As, re- 
spectively. The parameter values are taken as (a) gN{0) = 3, 
UN{0) = 0.5, (b) gN{0) = 3, UN{0) = 0.8, and (c) gN{0) = 2, 
!7iV(0) = 0.5. 



the s-wave state occurs. This result is consistent with 

the experimental fact that spin-triplet superconductiv- 
ity has not been observed in the ordinary metals. 

6.2 Narrow band systems 

Next, we examine systems in which the ratio Wc/tOc 
is not so large. From eq. (5.29), since Wc/u)c ^ 
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^/eF(W'^^7F)/u!B, such a situation is realized, for exam- 
ple, in semi-metals with extremely small cp and in narrow 
band systems with extremely small W. In such systems, 
since Wc ^ or Wc <C ep does not hold, the approxi- 
mation A''(e) w N{0) is not correct quantitatively, but it 
is not essential for our purpose. The phase diagram for 
q^p/2k^ = 0.5 and Wc/i^c = 5 is presented in Fig. 10. 
Since W > 2Wc = lOt^c ^ lO'^D, if we put uiu = 300 K, 
we obtain W > 0.3 eV. As the ratio Wq/uJc becomes 
smaller, since the retardation effect becomes weaker, s- 
wave pairing is suppressed more effectively by the on-site 
Coulomb repulsion. As a result, the area for the p-wave 
state is enlarged in the phase diagram, and Tcp becomes 
large enough to be observed in contrast to the result in 
the previous subsection §6.1. 

In Fig. 11, the transition temperatures for the param- 
eter values of Fig. 10 are plotted, where we have put 
(jJb = 300 K as an example. For UN{0) = 0.5, it is found 
that the transition temperature Tcp to the triplet state 
could reach 0.3 K at gN{0) w 2.6, while for gN{0) > 2.6, 
s-wave superconductivity occurs since T^s > Tcp. In con- 
trast, for UN{0) = 0.7, since Tcp > Tcs upto gN{0) w 4, 
because s-wave pairing is strongly suppressed by the on- 
site U. Thus, Tcp could reach a rather larger value 2 K, 
where the triplet state occurs. 

Figure 12 shows the dimensionless coupling constants 
Xa as functions of Wc/tOc- It is found that the p-wave 
state is favored for narrower band widths W 2Wc), 
because the on-site U suppresses s-wave pairing more ef- 
fectively. The coupling constants Ap for p-wave pairing 
depend on Wc/Wc slightly. If we put uJc ~ c^d = 300 K, 
we could estimate Tcp as Tcp = 0.015 K, 0.16 K, and 
0.79 K, for Ap = 0.1, 0.13, and 0.165, respectively. There- 
fore, for q^p/2kp = 0.5, which gives Atf ~ ^ 
o/tt ~ 0.3 X a, the transition temperature Tcp can be 
high enough to be observed in practice, where the s-wave 
state is suppressed by strong on-site Coulomb repulsion. 
For this mechanism, the smallness of the ratio Wc /u)c is 
essential. 

6.3 Systems with weak Coulomb screening 

Next, we examine systems with weak Coulomb screen- 
ing. As examined in our previous paper^^ and briefly 
discussed in §5.5, the screening can be weak, for example 
in layered systems with sufficiently large intervals be- 
tween the layers, although here we examine spherically 
symmetric systerms. 

Now, let us examine the case with q^p/2kp = 0.1 and 
Wc/ojc = 100. For q\-p/2k% = 0.1, we obtain Atf = 
V5/fcF > V5a/n - 0.7 X a. The short-range correlations 
lengthen the screening length by y/1 + nc as shown in 
eq. (3.18). In Fig. 13, the area of the p-wave state is much 
larger than those in Figs. 9 and 10. Thus, we confirm 
that as the screening becomes weaker, the anisotropic 
pairing interactions increase, as proposed in our previous 
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Fig. 13. The phase diagram for a weak screening ease, in which 
we set <j'^p/2fcp = 0.1 and Wq/oJc = 100. Along the vertical thin 
dotted line, the transition temperatures are plotted in Fig. 14. 




Fig. 14. The transition temperatures as functions of the effective 
constant g of the phonon-mediated interaction for a system with 
the weak screening, in which we set q^p/2kp = 0.1, Wc/u>c = 
100. We have put ujc = 300 K as an example. The solid and 
dashed lines show the transition temperatures Tcp and Tcs for 
the p-wave and s-wave states. 

In Fig. 14, the transition termpartures are estimated 

for the parameter values of Fig. 13. It is fomid that Tcp for 
the p-wave state could reach 0.3K, which are practically 
observable. 

6.4 Narrow band systems with weak screening 

Lastly, we examine systems with weak Coulomb 
screening and narrow electron bands. We set param- 
eters as q^p/2k^ = 0.1 and Wc/tOc = 20. The for- 
mar corresponds to the Thomas-Fermi screening length 
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Fig. 15. The phase diagram for a system with weak screening 
and a narrow electron band, in which we set (j';^p/2fcp = 0.1 and 
H'c/'^c = 20. Along the vertical thin dotted line, the transition 
temperatures are plotted in Fig. 16. 

Atf ^ 0.7 X a as mentioned in the previous case. The 
latter corresponds to VF/wd ~ 2Wc/wc = 40. If we 
put = 200 ^ 400 K as a typical value, wc obtain 
= 0.8 ~ 1.6 eV, which is realistic for the organic 
superconductors and the ruthenate superconductors. In 
Figs. 15 and 16, the phase diagram on the V-g plane 
and the transition temperatures are shown. Comparing 
Figs. 13 and 15, we confirm again that the area of the 
triplet state is enlarged when the electron band width 
becomes narrower. In Fig. 16, we find that the transi- 
tion temperature Tcp of the triplet state can reach 3 K 
at .gAf(O) w 6.2, ahhough for ,giV(0) > 6.2, the p-wave 
state is hidden behind the s-wave state (Tcs > Tcp). We 
shall discuss the application of this result to the real 
compounds in the last section. 

7. Summary and Discussion 

7. 1 Summary of the results 

Now, we shall summarize the results. In §2, we have de- 
rived an effective Hamiltonian which is written in terms 
of the short-range (on-site) Coulomb interactions U and 
the long-range Coulomb interactions Vq, from a basic 
Hamiltonian. On the basis of this effective Hamilton- 
ian, wc have examined the spin and charge fluctuations 
and screened Coulomb interactions (§3.1), considering 
the two-particle vertex part within an RPA extended to 
the model including both U and Vq. The effective inter- 
actions between electrons are written in the simple forms 
as shown in eqs. (3.11) and (3.12). We have derived an 
expression of the dielectric function and the screening 
length (§3.2). It was found that the screening length is 
lengthened by the on-site correlations due to U. 

In §3.3 and §3.4, we have examined the Coulomb 




Fig. 16. The transition temperatures as functions of the effective 
coupling constant g of the phonon-mediated interaction for the 
system with the weak screening and the narrow electron band, 
in which we set g^p/2fep = 0.1, Wq/uic = 20. We have put 
ojc = 300 K as an example. The solid and dashed lines show the 
transition temperatures Tcp and Tcs for the p-wave and s-wave 
states. 

screening corrections to the electron-phonon interactions 
and the phonon Green's function within the RPA con- 
sistent with those for the interactions between electrons. 
It was found that the sound velocity increases by the 
weakening of the screening due to the on-site correlations 
mentioned above. 

In §4, we have examined superconductivity. Wc have 
derived a general form of the effective Hamiltonian which 
includes the interactions mediated by the phonons and 
those mediated by the spin and charge fluctuations, 
where the short and long-range parts of the Coulomb 
interactions and the screening effects are taken into ac- 
count. It was found in eq. (4.15) that the effective inter- 
action is written in a summation of two terms, that is, 
and the renormalized phonon-mediated pairing in- 
teraction. Here, we call T^^ the Hubbard term, because 
it is nothing but the vertex part in the pure Hubbard 
model without the long-range Coulomb interaction Vq. 
In contrast, the latter term has the same form as the 
phonon-mediated interaction in the absence of the short- 
range Coulomb interaction U, but the prefactor and the 
phonon energy ^'q include the corrections due to U. 

In §5, we derived the effective interactions along the 
same way as the traditional weak coupling theory of su- 
perconductivity. In partictdar, for the system away from 
any magnetic instability, we derived the effective Hamil- 
tonian, with which we have examined anisotropic super- 
conductivity in §6. The present model includes the mod- 
els examined so far^"^' in appropriate limits except the 
corrections due to U. 

In §6, we have applied our theory to some ideal sys- 
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terns. The results are summarized as follows: (1) The 
triplet superconductivity occurs in the system with 
strong on-site correlations and weak electron-phonon 
coupling; (2) Spin-triplet pairing is favored when the 
band width is narrow and when Coulomb screening is 
weak;^° (3) The transition temperature of p-wave pair- 
ing is limited, because s-wave pairing is favored if the 
electron-phonon coupling is strong. In the most favorable 
system, the transition temperature T^p can be only of the 
order of 1 K as a crude estimation. These results agree 
with our knowledge as follows: The triplet superconduc- 
tivity has not been observed in the ordinary metals, while 
there are candidates in the heavy fermion superconduc- 
tors, and the layered superconductors with large layer 
intervals, such as (TMTSF)2X and Sr2Ru04, but their 
Tc's are not so high. 

7.2 Narrow band width and triplet pairing 

As mentioned in the above summary, we have found 
in Figs. 10 - 12 that p-wave pairing is favored in nar- 
row band systems. This result can be explained as fol- 
lows. The on-site Coulomb repulsion affects only s-wave 
pairing, but it has a limit itself. The effective Coulomb 
parameter /xq satisfies an inequality 

* ^ Mo < 1 ^ 

1 + /xo ln(Wc/a;c) ~ 1 + ln(Wc/a;c) ~ 

(7.1) 

where we have used 1 > UxoiQ,^) ^ UN{0), which is 
the necessary condition for the absence of the magnetic 
long-range order within the RPA. The upper bound /Umax 
is determined only by the ratio Wc/uic, which is usually 
of the order of W/uJu- Therefore, we obtain /imax ^ 0.18, 
0.30, and 0.38, for Wc/tOc ~ 100, 10, and 5, respec- 
tively. The phonon-mediated pairing interaction origi- 
nally has the anisotropic components, such as Ap, but 
they are much smaller than the isotropic component Ag. 
For triplet pairing to occur only by phonon-mediated 
pairing interaction, it is needed that ,s-wavc pairing is 
suppressed for some reason, and the dominancy of the 
coupling constants, Xg and Ap, changes. Prom the results 
shown in Fig. 12, the value /imax ^ 0.18 for Wc l^c ^ 100 
is too small to change the dominancy, when Atf ~ l/Ap. 
In contrast, the large value of /imax = 0.30 ~ 0.38 for 
Wc/loc = 10 ~ 5 can be large enough to suppress the 
s-wave state. 

7.3 Application to the heavy fermion systems 

The value of Tc obtained in Fig. 11 is consistent with 
Tc « 0.5 K in the heavy fermion UPts system,^^'^^ 
if the set of the parameter values q^-p/2k'^ = 0.5 and 
Wc/i^c = 5 is regarded to be appropriate for this com- 
pound. In addition, as in UGe2,^'^ if singlet pairing is 
suppressed by the internal exchange field due to the co- 
existing ferromagnetic long-range order, the p-wave tran- 
sition temperature T^p can be of the order of 1 2 K, 
which is consistent with the experimental value and our 



previous result.'' 

7.4 Conclusion and future study 

In conclusion, we have examined a model of the 
electron-phonon system with both the long-range and 
short-range Coulomb interactions, which is derived from 
a more basic model. In particular, we have examined su- 
perconductivity, and shown that the phonon-mediated 
interaction could induce spin-triplet pairing, when the 
Coulomb screening is relatively weak and the on-site 
Coulomb interaction is strong. The transition temper- 
ature can be high enough to be observed in practice, 
and similar to the experimental values in the candidates 
of the spin-triplet superconductors. In the future, more 
quantitative study on the isotope effect, in which the cor- 
relation effects and the strong coupling effects are taken 
into account, may give a direct evidence of the applica- 
bility of this mechanism. 

For more precise estimations of Tcp in the applica- 
tion to Sr2Ru04 compound, we need to take into ac- 
count the two-dimensionality. We have examined a two- 
dimensional effective model in our previous paper, and 
shown that a realistic value of T^p can be reproduced. It 
is possible to derive the effective Hamiltonian from the 
present Hamiltonian with the interaction of eq. (5.4). The 
derivation will be presented in a separate paper. "'^ 

The application to the quasi-one-dimensional (QID) 
organic superconductor (TMTSF)2X is another interest- 
ing subject to study as discussed in §1. In particular, 
Suginishi and the present author found by detailed cal- 
culations based on the present theory that specific fea- 
tures of the QID system play an important role in the 

pairing symmetry. It will also be presented in a separate 
97 

paper. ^' 
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